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ABSTRACT 

Recent observational studies have demonstrated that the majority of satellite galaxies 
tend to orbit their hosts on highly flattened, vast, possibly co-rotating planes. Two 
nearly parallel planes of satellites have been confirmed around the M31 galaxy and 
around the Centaurus A galaxy, while the Milky Way also sports a plane of satellites. 
It has been argued that such an alignment of satellites on vast planes is unexpected in 
the standard (ACDM) model of cosmology if not even in contradiction to its generic 
predictions. Guided by ACDM numerical simulations, which suggest that satellites 
are channeled towards hosts along the axis of the slowest collapse as dictated by the 
ambient velocity shear tensor, we re-examine the planes of local satellites systems 
within the framework of the local shear tensor derived from the Cosmicflows-2 
dataset. The analysis reveals that the Local Group and Centaurus A reside in a 
filament stretched by the Virgo cluster and compressed by the expansion of the Local 
Void. Four out of five thin planes of satellite galaxies are indeed closely aligned with 
the axis of compression induced by the Local Void. Being the less massive system, 
the moderate misalignment of the Milky Way’s satellite plane can likely be ascribed 
to its greater susceptibility to tidal torques, as suggested by numerical simulations. 
The alignment of satellite systems in the local universe with the ambient shear field 
is thus in general agreement with predictions of the ACDM model. 


Keywords: galaxies: haloes - formation - cosmology: theory - dark matter - large- 
scale structure of the Universe 


1 INTRODUCTION 


The Copernican principle dictates that the universe, at 
large, is homogeneous and isotropic. Yet the environment on 
galactic scales is anisotropic since galaxies reside within an 
complex network of filaments and sheets often dubbed “The 
Cosmic Web” ( [Bond et al.||1996| ). This anisotropic nature of 
the galaxy distribution in the Milky Way’s (MW) vicinity 


has long been evident since Kunkel & Demers (19761 and 
|Lynden-Bell| ( |1976[ |1982[ ) noted that the brightest satellites 
of the MW are distributed in a plane that is roughly per¬ 
pendicular to the MW disc (although there are difficulties 
in assessing this planar structure due to obscuration from 
the Galactic disc, e.g. |Willman et al.|2004[ ). The highly flat¬ 
tened plane of low-luminosity dwarf satellites seen around 
the MW was considered unique until these magnitudes could 
be probed elsewhere. [McConnachie et al.| ( |2009| ) recently 
compiled a photometric and spectroscopic survey known 
as Pan-Andromeda Archaeological Survey (PAndAS) of the 
stellar halo of M31, our nearest large galaxy, and uncovered 


a rich satellite system not unlike our own (see also Conn 
et al.|2011| 2012[). These observations revealed that roughly 


half (13 out of 27) of M31’s satellite galaxies belong to a 
“vast, thin, co-rotating” plane (Ibata et al.|2013 Conn et al. 


20131 viewed nearly edge on from the Earth. Tully (20131 


further suggested that a trio of dwarfs beyond M31’s virial 
radius (IC1613, ICIO and LGS 3) also aligned on the same 
plane found by Ibata et al. Follow up analysis by |Shaya fc| 
|Tully| ( |2013| of the remaining 14 satellites of M31 revealed 
that at least eight of them were also located on a second 
plane (although not necessarily a co-rotating one). The two 
planes of satellites around M31 are roughly parallel to each 
other, offset by around 15 degrees. 


A bit further afield, Centaurus A (hereafter CenA), a 
massive galaxy with around 30 dwarf satellites, shows a very 
similar set up to M31: [Tully et al)| ( |2015[ ) found that all but 
2 out of 29 of the CenA satellites, for which distances are 
known, appear to cluster onto two roughly parallel planes 
(offset by 8 degrees) each incredibly thin (r.m.s. values of 
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~50kpc), fairly long (extending ~ 250kpc) and separated by 
around 300 kpc. Beyond these three systems, satellite galax¬ 
ies have been found to exhibit flattened distributions around 
M81 ( |Chiboucas et al.|2013| ) and NGC 3109 ( [Bellazzini et al.| 


20131. Using stacking techniques, similar anisotropic distri¬ 


butions of satellites have been found in large surveys (Zarit- 
|sky et S^|1997l such as the 2dF (Sales & Lambas 20041 


and SDSS (Brainerd 2005 Yang et al. 2006 Agustsson & 
Brainerd|2010 (. 


Although not identically defined, the evidence for the 
ubiquitous anisotropic distribution of satellite galaxies has 
by now a strong observational foundation. The abundance 
of evidence pointing towards the anisotropic distribution of 
satellites around large central host galaxies calls for a com¬ 
prehensive and robust explanation, within the framework 
of the cold dark matter with dark energy (ACDM) model. 
Accordingly gravitational instability drives the formation of 
structure via the amplihcation of small primordial fluctua¬ 
tions that prevailed in the infant universe ( Zerdovich|1970 1. 
In doing so, small initial deviations from perfect spherical 
symmetry are also amplihed. The cosmic web is a manifes¬ 
tation of the anisotropic nature of the instability. 

ACDM makes very strong predictions on how the large- 
scale structure evolves in general, and fully accounts for the 
collapse of dark matter halos. However its constraints on 
the distribution of baryons, in particular of stars and stellar 
systems, are considerably less conclusive. On small scales, 
this is exemplified by the so-called “missing satellite” prob¬ 
lem ( [Klypin et al.[|1999 Moore et al.|1999 1 and the “too big 
to fail” problem Boylan-Kolchin et al. (20111. A number of 
solutions to both of these issues have been postulated, many 
of which rely on correctly approximating the complex inter¬ 


play between baryons and dark matter (e.g. Di Cintio et al. 
|2011[ [Sawala et al.|[^14[ ). It is by now clear, that baryons 
must play a crucial role on these scales. 

Accordingly, gas-dynamical numerical simulations (e.g. 
|Schaye et al. |2015[) and A^-b ody runs with semi-analytical 
models (e.g. Springel et al. 20051 have achieved unprece¬ 
dented success in describing the galaxy population as a 
whole. However such models have continued to struggle to 
explain the origin of satellite galaxy setups such as those seen 
in the local universe (although have not altogether failed, 
see |Kang et al.||2005| [Zentner et al.||20d^ [Libeskind et ah] 


2005 


2007 


Deason et al. 20111. Since most ACDM simu¬ 


lations have consistently found it difficult to explain these 
satellite planes ( jPawlowski et aL]|2012[ |2014| , a number of 
authors have suggested alternative formation mechanisms, 
namely that they are of tidal origin (e.g. Kroupa et al.|2005{ 


jPawlowski et al^|2011[ jMetz fc Kroupal 


20071, possibly en¬ 


gendered by ancient merger with M31 ( Fouquet et al.|2012 |. 

A less controversial issue is how and to what extent 
the cosmic web imprints its preferred directions on galax¬ 
ies in general (e.g. jAubert et ^|2004| ). Recently, jLibeskindj 
jet al.j ( |2014b| ) showed that satellites are channeled towards 
their hosts along directions dictated by the large scale struc¬ 
ture, irrespective of environment, redshift, or host/subhalo 
mass. Specifically the eigenvectors of the velocity shear ten¬ 
sor, used also in this work, determine to a large extent, the 
orientation of the preferred axes along which satellites are 
accreted (see also jLibeskind et al.|20lT| . 

Motivated by these studies which link the large and 
small scales of structure formation, we aim to examine the 
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Table 1. Directions in the Local Volume in supergalactic coordi¬ 
nates. The upper rows include, rny and fyirgo- the unit vectors 
in the directions of the Local Void and the Virgo cluster, iimsipi 
and ri]y[ 3 iP 2 ' unit normals to the two planes of satellites around 
M31. hcAPl and hcAP 2 : unit normals to the two planes around 
CenA. hjviwP- unit normal to the MW plane. The bottom rows 
include each of the three components of the orthonormal eigen¬ 
vectors (ei, 02 , and 03 ) of the shear field are evaluated at the 
Local Group and at the location of Cen A. Note that the eigen¬ 
vectors are non-directional lines and the +/- direction is thus ar¬ 
bitrary. References for the directions in the top rows are: [ 1 ][T^ 
et al. IpOO^; [2] NED database http://ned.ipac.caltech.edu/, [3] 
Shaya &: Tully|2013| Table 2; [4] [Tully fc et|201^ [5] jPawlowskij 
Hz Kroupa|2013| 


relationship between planes of satellites and the velocity 
field in the local nniverse. The distribntion of the satellites 
around host galaxies in general, and around the MW, M31 
and CenA in particular, raises one main issues: can such 
thin partially co-rotating setups be accommodated within 
the ACDM model or are these anisotropic distributions at 
odds with generic predictions? Is the observed large-scale 
structure velocity field responsible for the peculiar satellite 
alignments seen around the MW, M31 and Cen A? 


2 METHOD 


Numerous schemes for characterizing the galaxy distribu¬ 
tion on the Megaparsec scale have been suggested, each with 
their own pros and cons depending on the task at hand. In 
this work, we describe the cosmic web using peculiar veloc¬ 
ities (the dynamical contribution due to the general expan¬ 
sion is not considered). In order to do so we appeal to the 
Cosmicflows-2 (CF2, jTully et al. 20131 survey of peculiar 
velocities. Presently this is the richest and deepest full sky 
survey of its kind. It enables the reconstruction of the under¬ 
lying density and 3D velocity fields out to distance exceeding 
150 Mpc. The quality of the data allows for the mapping of 
the nearby universe in great detail: in particular it has es¬ 
tablished the existence of the Laniakea Supercluster that 
includes the Local Croup and the entire Local Supercluster 
( jTully et aLTMIl ). 
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2.1 Wiener Filter reconstruction 


The 3D velocity field and its associated shear tensor are 
calculated from the Cosmicflows-2 dataset by means of 
the Wiener Filter (WF) and constrained realizations (CRs 
Zaroubi et al.|ll)OT Courtois et al.|2012 20131. The ACDM 
model with WMAP5 ( [Komatsu et al.| 2009| cosmological 
parameters (i.e. ilo = 0.228, Q,b = 0.046, Ha = 1 — flo, 
Ho = lOOh = 70 km/s/Mpc, as = 0.812) is used as the 
prior model for the reconstruction. The reconstructed WF 
density and 3D gravitational velocity fields are evaluated in 
a box of 320 Mpc//i side-length. A Fast Fourier Transform 
(FFT) is used to construct the CRs, hence periodic bound¬ 
ary conditions are assumed. The construction of the CRs 
has been done with a box of 640 Mpc//i side-length, to sup¬ 
press periodic boundary conditions artifacts. Both WF and 
CRs fields are spanned on a grid of 256®. The WF provides 
the mean field given the data and the prior model (the CF2 
dataset and the ACDM/WMAP5 model in our case). 

In order to estimate the robustness (or accuracy) of the 
WF reconstruction, an ensemble of 20 CRs is then used to 
sample the scatter around the mean field. In particular it 
is used to evaluate the error-bars of the eigenvalues (ie not 
their magnitudes) and scatter in direction of the eigenvec¬ 
tors of the shear tensor. We note that the spatial resolution 
of the CRs is double that of the WF, and therefore the un¬ 
certainties derived from the CRs are evaluated on the scale 
of 5.0 Mpc/h, compared with the 2.5 Mpc/h of the WF. As 
such we do not expect the magnitude of the eigenvalues of 
the two methods to be in agreement, rather the method is 
used to get a handle on the error bars. 


2.2 Shear Tensor 


Spatial derivatives of the reconstructed peculiar velocity 
field are taken and the dimensionless shear tensor is defined 

as Eq /3 = where a and l3 are the x,y, 

and 2 components of the velocity v and position r and Ho 
is Hubble’s constant. A minus sign is added for convention. 
The tensor’s sorted eigenvalues (Ai > A 2 > A 3 ) describe 
the strength of compression (positive values) or expansion 
(negative values) along the eigenvectors of the shear: ei, 
e 2 and 63 (see Hoffman et al. 2012[ Libeskind et al. 2013 


|2014a| ).The spatial derivatives that are used to build the 
shear tensor are evaluated by the FFT, hence the formal 
resolution of the shear is 2.5 Mpc/h for the WF and 5.0 
Mpc/h for the CRs. It follows that in the present work the 
entire Local Group (i.e. the MW and the M31) experience 
the same velocity shear, nearly identical to the one experi¬ 
enced by the adjacent Cen A (see below). 

The eigenvector directions that we accept in this work 
are compared with published values for the normals to the 
five satellite planes found around the MW, M31 and Cen 
A as well as other objects in the local universe. These are 
presented in the top section of Table. 


3 RESULTS 

A Wiener Filter reconstruction applied to the CF2 data re¬ 
covers the underlying 3D gravitationally induced velocity 
field with an effective resolution of a few Megaparsecs. At 
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Table 2. Alignment between the shear tensor and relevant direc¬ 
tions in the local volume. We examine primarily ei, the axis of 
fastest collapse, but also show how the local filament, defined by 
03 , aligns well with the direction of the Virgo Cluster. fi,v and 
fvirgo are the unit vectors in the directions of the Local Void and 
the Virgo cluster, umsipi and fiM 3 iP 2 are the normals to the two 
planes of satellites around M31. Similarly fiQAPl and ncAP 2 are 
the normals to the two planes around CenA. Amwp is the normal 
to the MW satellite plane. 


the location of the Local Group and Gen A, the eigenval¬ 
ues are found to be (Ai, A 2 , A 3 ) = (0.032, 0.018, -0.020) 
and (0.029, 0.011, -0.025), respectively. Considerable infor¬ 
mation is encoded in these eigenvalues and their relative 
magnitudes. The sign, combined with the fact that A 2 is 
much closer to Ai than it is to A 3 suggests that the Lo¬ 
cal Group and Cen A reside in a filament compressed along 
two eigenvectors (ei and 02 ) and expanding along the third 
(es). This picture is also evident in the spatial cosmography 
and reconstructed velocity field, shown in Fig. The shear 
eigenframe at the location of the Local Group is indicated by 
the three arrows (the 62-03 plane is shown by a grey disc). 
The shear eigenframe does not change much 3.8 Mpc away 
at the location of Cen A. It is clear from both the velocity 
field as well as the density field that both systems reside in 
a filament, with Cen A being right on the border between 
a small void and the local filament. The stretching induced 
by Virgo along 63 and the compression due to a void in the 
direction of 62 are clearly manifested. The Local Void is here 
above the slab, in the positive supergalactic z direction, and 
its compression is directed along ei. In the bottom section 
of Table we present the eigenvectors (ei, 62 , 63 ) of the 
velocity shear tensor experienced by the Local Group and 
Cen A, in supergalactic coordinates. 

The shear tensor eigenvectors provide the natural prin¬ 
cipal frame of reference within which the anisotropies and 
preferred directions of the environments of the Local Group 
and Cen A can be studied. How these planes align within 
the shear eigen-frame, can be quantified by examining the 
angle made by each plane’s normal and the eigenvectors. 
Fig. [2] shows how the five planes of satellites in the local 
universe (the MW satellite plane, and the two around M31 
and Cen A) are oriented with respect to the principal di¬ 
rections defined by the local shear tensor. The relationships 
are quantified in Table. and described here. The direc¬ 
tion towards the center of the Local Void is aligned with 
61 . The spine of the local filament ( 63 ) points close to the 
direction of the Virgo Cluster. Most interestingly, the local 
flow dictates the preferred directions of the M31 and Cen A 
satellites. The normals to two planes around M31 and Cen 
A lie very close to 61 , all within ~ 15 degrees (with one of 
Cen A’s planes being within ~ 1 degree!) The normal to 
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Figure 1. The cosmic peculiar velocity and related over-density field in the local universe, reconstructed from the CosmicFlows-2 survey, 
by means of the Wiener Filter methodology. In both plots stream-lines are colored according to the local speed, shown by the color bar 
at right. In the top plot (a), the over-density field is colored from the highest density regions (the Virgo Cluster in red) to under-dense 
ones (voids, in blue). In the bottom plot (b) all galaxies from the V 8 k survey | |Tully et al.||200^ in this 20 x 20 x 3Mpc//i slab are 
plotted as black dots. The flow towards the Virgo Cluster, as well as the “local filament” that bridges the Local Group with the Virgo 
Cluster, exhibits clear features in both the density and velocity fields. Centaurus A is located on the edge of this filament. The Local 
Void, which sits above the slab in the positive SGZ direction, pushes down onto the local filament. The shear eigenframe centered on 
the LG is indicated by the three eigenvectors. The ( 02 - 63 ) plane is shown as the grey disc in (b). The ei eigenvector points mostly out 
of the slab. Note how the local filament and flow toward Virgo are aligned closely with 63 . In both figures, labels for some well known 
objects are also included. 


the MW satellite plane is not as well aligned and is offset 
by around 38 degrees. We provide a possible explanation 
for this mis-alignment in section | 3 . 2 [ The close alignment 
of four out of five satellite planes of the local universe with 
the (02-03) plane is clearly visible. It should be added that 


the line connecting the MW to M 31 lies within the (02-63) 
plane, well away from 01 (and ~40 degrees away from 03). 

The following picture emerges: the Local Group and 
Cen A reside in a filament defined primarily by the dark 
matter distribution. This filament is stretched by the Virgo 
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• ^cAPi Cen A plane 1 
^CAP 2 tp Cen A plane 2 

* n„wp to MW sat plane 

^ ^Locol Void 
B ^Vlrgo 

— M31 Plane 1 

— M31 Plane 2 

— Cen A Plane 1 
Cen A Plane 2 
MW satellite plane 





Figure 2. Satellite planes in the local volume projected in the eigenframe of the shear field show a clear tendency to align with the 
62-^3 plane as shown here in an Aitoff-Hammer projection of the unit sphere. Colored curves denote the planes themselves with a pair of 
corresponding symbols denoting the it directions of the normals to each plane. The ei, e 2 , and 03 axes are indicated by the labels. Black 
curves around these labels show regions of 15 degrees. The angles formed between these satellite planes and the shear eigenframe are 
given in Tab.[^ Four out of five satellite planes in the Local Volume are well aligned with their local shear field. Note that one Centaurus 
A satellite plane (CAP2, cyan curve) is within less than a degree of the 02-63 plane. The Milky Way satellite plane on the other hand, is 
rotated about the 62 axis by around 38 degrees. This offset can be attributed to a smaller Milky Way halo mass. The direction toward 
the Local Void (cLocal Void) is also close to ei, while the direction toward Virgo (rvirgo)i viewed either from Cen A or the Local Group, 
is close to 63 . 


Cluster, while being compressed from the supergalactic pole 
by the Local Void and from the sides by smaller voids. The 
four satellite planes that surround the two most massive 
galaxies, M31 and Cen A, lie within or close to the ( 02 - 
03 ) plane. Both of the M31 satellite planes, the two Cen 
A satellite planes, and the MW-M31 orbital plane are, to 
within observational uncertainties for the Local Group spin, 
roughly parallel with each other and with the ( 02 - 63 ) plane. 
It follows that influences on scales of tens of Megaparsecs, 
notably the Virgo Cluster and the Local Void leave their 
imprint on the anisotropies that characterize the satellite 
distribution around Cen A and M31. 


3.1 Statistical significance in the estimated shear 
eigenvalues and eigenvectors 

As mentioned in section [27T| we estimate how robust the de¬ 
termination of the WF shear field is by using CRs. In the 
standard model of cosmology the over-density and velocity 
fields, on large enough scales, constitute random Gaussian 
fields. It follows that the individual matrix elements of the 
shear tensor constitute random Gaussian variables whose 
mean value is zero. The eigenvalues of the tensor are not 
Gaussian variables, and therefore the mean value of each 
eigenvalue of the ensemble of GRs differs from the corre¬ 
sponding eigenvalue calculated from the WF reconstruction. 
In particular the mean amplitude of the GRs eigenvalues is 
expected to be larger than for the WF. The WF provides 
a more conservative estimation of the eigenvalues than that 
of the GRs and the WF value is used here to characterize 
the environment of the Local Group. The GRs are used to 
estimate the statistical uncertainties for the WF, which is 
presented in section [XT] 

The shear tensor has been evaluated for the 20 GRs. 
Table presents the mean and the standard deviation of 


quantity 

mean value ± std. dev. 

Ai 

0.148±0.038 

A 2 

0.051±0.039 

A 3 

-0.160±0.033 

Ai — A 2 

0.097±0.042 

CO 

1 

0.211±0.061 

eWF.eCR 

0.890±0.137 

pWF . pCR 

0.892±0.136 

WF . CR 
*^3 

0.923±0.047 


Table 3. Stability of the Wiener Filter reconstruction: The mean 
value ± the standard deviation of the eigenvalues, the eigenvalue 
difference and the alignment of WF and CR eigenvectors. Given 
the different smoothings in the two methods (CR versus WF) we 
do not expect the eigenvalues to have the same magnitude rather, 
its the small value of the standard deviation which is important. 

the eigenvalues (Ai , A 2 , A 3 ), as well as the differences in 
eigenvalues (Ai — A 2 and A 2 — A 3 ) and the stability of each 
eigenvector (e)'^^ • ep^). The filamentary nature of the lo¬ 
cal environment is thus confirmed at around the 3a level. 
The table further shows the mean and standard deviation 
of the alignment of each eigenvector of the GRs with the 
corresponding eigenvector of the WF field. The stability in 
the direction of each eigenvector is statistically very signifi¬ 
cant. Note that the direction of the filament, given by 63 , is 
statistically the most robust eigenvector. 


3.2 The case for the MW mis-alignment 

As shown in Fig. the MW satellite plane is not as well 
aligned with the shear field as these other four systems. The 
normal to this plane appears to have been rotated (around 
62 ) by about 38 degrees. This offset may not be a surprise if 
MW is less massive than its partner M31 as suggested by a 
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number of studies (Smith et al.|20d7 Xue et al.|2008 |Kafle| 


|et al.|2014| |Penarrubia et al.|20i4 Piffl et al.|2014a|b| ) This 

is because the satellite system of a less massive halo in a 
pair like the LG, is dynamically less stable and more prone 
to torques, in this case that act in the ei-ea plane. 

In order to gauge how the alignment of a plane of satel¬ 
lites with the shear eigenframe may depend on the host halo 
mass, we turn to numerical simulations, specifically the Bol¬ 


shoi simulation (Klypin et al. 20111, a large cosmological 


box (of periodic side length L = 250Mpc//i) populated with 
2048® dark matter (DM) particles, each with a particle mass 
of ~ 1.9 X 10® Mq. The simulation was run with cosmological 
parameters that are consistent with WMAP5 and the 2 = 0 
snapshot is considered. After particles are grouped into DM 
haloes by use of the BDM halo findeij^ ([Klypin &: Holtzman 


1997 Riebe et al.||2013 l, we divide the halos into two sub¬ 


samples that mimic pairs of Local Group galaxies and Gen 

A. 

To identify Local Groups, we search for all pairs of 
galaxies: (1) with masses between 5 x 10^® and 5 x 10®®Mq; 
(2) with a separation between 0.5 Mpc and 1.5 Mpc; (3) 
that have no halo with mass > IO^^Mq within 10 Mpc of 
the pair barycenter; (4) and that have no other halo with 
mass > IO^^Mq within 0.5 Mpc of either candidate. To iden¬ 
tify Gen A look-a-likes we simply find all haloes (1) whose 
mass is between 5 x and 5 x IO^^Mq and (2) that 

have no mass greater than itself within 3 Mpc. These criteria 
are designed to be fairly general so as not to impose strong 
a priori constraints on any subsample, yet are designed to 
match some of the more accepted characteristics of these 
galaxies. 

Once the “Cen A” and “Local Group” samples are 
amassed, their shapes are computed as defined by the iner¬ 
tia tensor. The shortest axis of the ellipsoid Chaio is assumed 
to point in the same direction as the normal to the plane of 
satellites. The shear tensor is then computed from the ve¬ 
locity field on a regular grid with a 2.5Mpc smoothing (see 
previous work, Libeskind et al.|2013 1 and the angle between 
Chaio and the eigenframe is examined. In Fig. [^ we present 
the distribution of the cosine of the angle formed between 
Chaio and © 1 . 

Fig.jSj shows that the alignment between the less mas¬ 
sive Local Group partner (red), the more massive Local 
Group partner (black) and isolated Gen A-like halos (blue). 
The dotted curves show the alignment of the two Local 
Group partners when only pairs with a mass ratio of less 
than 1:2 are considered (motivated by jPefiarrubia et J)] 


20141. A number of interesting points can be inferred here. 


Isolated massive Gen A-like haloes are well aligned with the 
shear field - their Chaio axes are roughly 25% more likely to be 
tilted towards the ei axes than expected from a random dis¬ 
tribution, making them inconsistent with a uniform distribu¬ 
tion at the 6 ct level. Local Group partners, too, show a ten¬ 
dency to have their Chaio axes aligned with the ei axis. The 
more massive Local Group partner (M31-proxy) is around 
18% more likely than a random distribution to be aligned 



Icos 01= . e, 


Figure 3. The alignment of dark matter halos with the shear field 
in a large cosmological simulation shows similar behavior to that 
observed in the local universe. The probability distribution func¬ 
tion of the cosine of the angle formed between the short axis of a 
halo (chaio) and the axis of fastest collapse (ei) of the surround¬ 
ing shear field is shown for isolated massive Cen A-like haloes 
(blue) and for the more massive (black) and less massive (red) 
partners in LG-like set ups. The dotted line shows the alignment 
of each LG partner when only LGs with mass ratios less than 
1:2 are considered. The shaded gray region represents the scatter 
about = 1 that one would expect from a 100,000 uniform 

distributions of the same size (~50,000 pairs), thus making these 
probability distributions highly incompatible with random ones. 


with the shear eigenframe. The less massive partner (MW 
proxy) has a significantly weaker alignment with the shear 
eigenframe, with just 14% better aligned than a random 
distribution. This result is true for all mass ratios, but is ex¬ 
acerbated for Local Groups where the mass ratio is less then 
1:2 (dotted curve). Indeed, with these mass ratios, the more 
massive partner is around 20% more aligned than a random 
distribution; the less massive partner is just 9% more aligned 
than a uniform distribution. 

The conclusions we arrive at from looking at how DM 
halo shapes are oriented with respect to the shear field are 
clear: the more massive the object, the better the alignment; 
the more isolated the halo, the better the alignment; in Lo¬ 
cal Group like pairs, the more massive halo is better aligned 
than the less massive one; when the mass ratio is far from 
1:1, the difference in alignment for the more and less mas¬ 
sive partner grows. Indeed for mass ratios of 1:10 or more 
(not shown) the less massive Local Group partner shows no 
alignment at all. 


^ Much post-processed data from the Bolshoi (and other 
simulations) including halo and subhalo catalogues. Merger 
Trees, and velocity and density fields is publicly available: 
http: / / www.cosmosim.org/cms/simulations/multidark-project / 


4 SUMMARY AND CONCLUSIONS 

Vast planes of dwarf galaxies have been discovered through¬ 
out the local universe. First spotted around the MW, 
( jKunkel fc Demers|1976||Lynden-Bell||1976||1982||Pawlowski| 
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fc Kroupa|2013 1, a pair were found around M31 (|Ibata et al.| 


2013 Conn et al.||201^ [Shaya fc Tully|[2013|) and mo st re¬ 


cently, a pair was found near Centaurus A (Tully 2015 Tully 
|et al.|20T^. In the case of the MW and at least one of M31’s 


planes, the dwarfs appear to be co-rotating (Metz et al. 2008 
Pawlowski et al.|2014l. 


We have applied the Wiener Filter algorithm ([Zaroubi] 


et al.||l999 l to the CosmicFlows-2 survey of peculiar veloc- 
ities ( [Tully et ah ]| 20 T^ , to reconstruct the full density and 
velocity field in the local universe. Such a reconstruction has 
allowed us to examine how the five known planes of dwarfs 
are embedded in the cosmic web. We have found that the 
main cosmographic feature of the local universe is a filament 
that stretches from the Virgo Cluster, some 16 Mpc away, 
across the Local Group that is not seen in the galaxy dis¬ 
tribution but clearly evident in the velocity field and dark 
matter density field. Our main result is that four out of five 
of the planes of dwarf galaxies that have been discovered are 
well aligned with the main direction of compression, defined 
by the principal axis of collapse (ei) of the shear tensor. 
This direction is close to the direction of expansion of the 
local void. One of the dwarf planes, that of the Milky Way, 
is not as well aligned with the shear. In Section [3.2| we sug¬ 
gest that this misalignment may be attributed to its lower 
halo mass, and consequent susceptibility to torques which 
appear to have rotated the Milky Way satellite plane about 
the 02 axis by around 38 degrees. The main conclusions that 
we can thus draw, is that structure on the scale of ~ lOMpc 
appears to be having a direct effect on small scale structures 
that are on the order of the virial radius. 

How does such a cosmographic picture fit into the 
ACDM model? |Libeskind et al.| ( |2014b[ ) showed that in nu¬ 
merical simulations, the infall direction of sub-halos onto 
their host halos is determined by the eigenvectors of the ve¬ 
locity shear tensor (ei, 02 , and 03 ). Subhalo accretion dis¬ 
plays a clear tendency to occur well away from ei, the axis of 
fastest collapse, and preferentially in the ( 02 - 63 ) plane with 
a propensity to be aligned with 63 . This accretion pattern 
extends over all redshifts and over scales much larger than 
the virial radius of a given halo. Furthermore it happens 
irrespective of the host or subhalo mass, or environment: 
namely in voids, sheets, filaments, or clusters (hence it is 
known as “universal accretion”) - it is characterized only by 
the directions of the shear eigenvectors. Alignments of sub¬ 
haloes with the cosmic web in the Sloan Digital Sky Survey 


tend to support the idea of universal accretion (Ibata et al. 
|2014[|Lee fc Choi|2015j |Tempel et al.|2015[ ). 

Although weakened by processes associated with virial- 
ization, the consequence of universal accretion is also seen in 
the present epoch distribution of sub-halos within their hosts 


in numerical simulations such as Libeskind et al. (20121 or 
|Shi et al.lpoTs] ). Such flattened distributions as seen in the 
real universe are not readily found but are to some extent 
reproducible if the criteria to define them is kept vague and 


not bound by so-called a posteriori statistics (e.g 

Libeskind 

et al.|2005| Zentner et al.||2005| Libeskind et al.|2009| 

Dea- 

son et al. |2011 Sawala et al.|2014 Lovell et al.| 

2011 

<ang 

et al.|2005| Wang et al.|20131 [Cautun et al.|2014| 

Cillet et al. 

20151. A number of these studies, admittedly, find that the 


criteria needed for a successful z = 0 match are likely to only 
produce transient, unstable planes when examined at higher 
2 . What is much more stable is the alignment of haloes with 


the eigenvectors of the cosmic web. Spin and halo shape are 
directly built by the stresses and strains due to the grav¬ 
itational deformation of the cosmic velocity field. This is 
clearly seen in numerical simulations of structure formation 
and now for the first time demonstrated for the local uni¬ 
verse, a region whose velocity field we are able to probe with 
high accuracy. 

That satellites occupy primarily the ( 62 - 63 ) plane can 
be thought of as a general consequence of the formation 
of “Zeldovich pancakes” (Zel’dovich 19701, wherein matter 
first collapses into vast cosmic sheets, defined primarily by 
61 . On smaller scales, these sheets are composed of small 
clumps that hierarchically merge according to the “bottom- 
up” scenario, suggested by the ACDM model ([White fc Rees| 


19781. Since initial collapse has already occurred along 61 , 


not much material is left to be accreted from this direction. 
Instead haloes are guided along the ( 62 - 63 ) plane, towards 
density peaks where their ultimate fate is to merge. 

While an examination of the directions in the local uni¬ 
verse has led us to an understanding of how satellite planes 
are aligned, a number of open issues remain. What deter¬ 
mines the extreme thinness and/or apparent co-rotation of 
some of these planes? Are the planes seen in ACDM simula¬ 
tions stable or do they precess and evaporate in time? Also, 
the azimuthal distribution of M31’s satellites is lop-sided, 
with many galaxies found on the near side, closer to the 
MW. Are these galaxies reacting the the MW’s presence, or 
does this aspect of the problem have a different origin? We 
leave these important challenges to future work. 
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